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Abstract. The paper develops a method for discrete computational Fourier analysis of functions defined 
on quasicrystals and other almost periodic sets. A key point is to build the analysis around the emerging 
QQ ' theory of quasicrystals and diffraction in the setting on local hulls and dynamical systems. Numerically 

f™*^ | computed approximations arising in this way are built out of the Fourier module of the quasicrystal in 

, question, and approximate their target functions uniformly on the entire infinite space. 

^s^J ■ The methods are entirely group theoretical, being based on finite groups and their duals, and they are 

practical and computable. Examples of functions based on the standard Fibonacci quasicrystal serve to 
OJJ 1 illustrate the method (which is applicable to all quasicrystals modeled on the cut and project formalism). 
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1. Introduction 

In this paper we consider the problem of discrete methods for dealing with functions that are intrinsically 
'-P' ■ almost periodic, but not actually periodic. Quasicrystals, quasicrystalline photonic crystals, Faraday wave 

experiments, and other physical phenomena arising from the interaction of incommensurate frequencies, all 
_P ■ display the features of almost periodicity. As a typical example one may think of a potential field of a physical 

quasicrystal. The salient features of quasicrystals are highly structured long-range order (represented by pure 
point or near pure point diffraction) but no periodic order. Thus the potential is not a periodic function, 
but rather belongs to the domain of almost periodic functions. 

Here we put forth a method for finite discrete analysis of almost periodic functions that has the following 
main features: 

■^j- ' • It is entirely based on group theoretical methods, primarily finite groups and their duals. 

• The discretely computed Fourier approximants are themselves almost periodic and uniformly ap- 
proximate their target functions over their entire domains. 

• The Fourier frequencies involved in the approximation lie in the module of Fourier frequencies of the 
target function. 

A standard approach to modeling such structures is to take a finite part of it, impose periodic boundary 
conditions, rationalize and reduce the object to a periodic approximant, and then apply usual crystallography. 
Although this type of periodization is used routinely and successfully for many modelling problems in the 
theory of quasicrystals, it is not entirely satisfactory. Almost periodic order goes beyond periodic order in 
fundamental ways, its essence appearing as a underlying incommensurability which pervades every part of 
the theory. For instance, a key feature of quasicrystals appears in the Fourier module which parameterizes 
the Bragg spectrum and always has rank higher than (typically double) the dimension of the ambient space 
of the quasicrystal. Periodization destroys this, and by its nature can only produce results that can fit data 
over the finite range specified by the imposed periodization boundaries, whereas the essence of the material 
is that its order is long-range. The present paper does not involve any periodization and avoids these issues. 

The theory of almost periodic functions was initiated by H. Bohr [7\, on the basis of earlier work on 
uniform approximation of functions by trigonometric polynomials by P. Bohl. It was greatly extended by 
the work of A. Besicovitch 3 ], S. Bochner, J. von Neumann, [KHE], N. Wiener [SJ, H. Weyl [53], and 
others [9j[15j[I]- The advent of quasicrystals and aperiodic tilings instigated a revival of the field, and led 
to extensive study of the cut and project formalism and the theory of pure point diffraction [16j H3 HU E2] 
which have become the mainstays of experimentalists and theorists alike. An important component of this 
is the use of dynamical systems and dynamical hulls. These are ideally suited to the phenomenon of almost 
periodicity, which appears in the dynamics as recurrence, and provide a natural setting for the Fourier 
analysis that is used to study it. 

Although our study is of almost periodic functions, their importance in the subject of quasicrystals is 
that they arise from functions whose behaviour is dominated by the local environments of the quasicrystal in 
question. The way in which this happens is mostly taken for granted, but in fact there are some interesting 
assumptions involved, and for this reason we begin by formally defining local functions with respect to a 
given point set A, and showing how it is that they are connected with almost periodic functions. 

l 
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We can create a dynamical hull (X, /x) from A. This is a compact space that arises from A and its 
translations, and it lies at the heart of the Fourier analysis of local functions on A. One assumes (in many 
important cases it is forced) a probability measure fj, on which M. d acts in a measure preserving way. A local 
function / lifts to some new function F on (X, /i), and this is an L 2 function. Now the analysis of / can be 
related directly to the analysis of F, and for this we have a powerful tool in the form of the action of M. d on 
L 2 (X, /i), which is unitary. All of this material is explained in §2. 

To go further, we next place ourselves in the situation of the cut and project formalism, which is the 
standard method of modeling used in the study of quasicrystals. The set A is now assumed to be a model 
set (cut and project set). In this setting the hull (X, ^), and more particularly L 2 (X,[i), can be described 
explicitly in terms of a higher dimensional torusQ(higher dimensional periodicity!), and it is a straightforward 
matter to carry out Fourier analysis of F. It is the restriction of this Fourier analysis back to / that provides 
the required Fourier analysis of /. By its very nature this is almost periodic and captures the full aperiodic 
nature of A, including the correct Fourier module in which physical information actually appears. This is 
the content of §3. Readers familiar with the cut and project method who do not wish to go through the 
theory of local hulls and local functions may read H3.ll and M3.4L and then continue starting at SjD 

This theoretical analysis is based on higher dimensional structures that are not explicitly computatable, 
as well as the usual array of countably many Fourier coefficients, each of which is the outcome of integration. 
To be a practical tool, the analysis has to be reduced to finitely many objects that are explicitly computable 
entirely in the context of the given function /. The resulting approximants are trigonometric polynomials 
(quasi-periodic functions) whose frequencies come from the Fourier module of the original function. In §4 
we outline the method of discretization, which depends primarily on the construction of a refinement lattice 
of the lattice of the cut and project scheme and, along with it, its dual lattice. Together these produce two 
finite groups which are in Z-duality to each other. The selection of data points and appropriate Fourier 
frequencies is governed by these two groups. 

This is best illustrated by examples and for this we have chosen two local functions based on one of 
the famous Fibonacci point sets. This has the advantage of being straightforward to construct and easy to 
visualize, while at the same time containing all the essential features of more general model sets. §5 prepares 
the mathematics of the Fibonacci cut and project scheme, and §6 shows some explicit computations for the 
particular local functions we have chosen. As is evident from the examples, the main effort required is in 
the creation of the data points. Once this is done, the same set of data points and Fourier frequencies will 
work for any almost periodic function arising from the same cut and project scheme. 

The results are striking in two ways. First of all, the approximating functions are remarkably good, 
given the amount of data from which they are produced. Second, the approximating functions are not just 
local approximations they are also global approximations, in the sense that they provide finite Fourier series 
that approximate / throughout its entire domain (namely R). Of course this was to be expected, but it is 
impressive to see it in action. 

The numerical methods we introduce here are designed to be efficient and, of course, to utilize the 
inherent almost periodicity. In the development of the theory we provide error estimates that can be worked 
out specifically in the cases of interest. We note that the primary weakness in the error estimates is not one 
that is due to the aperiodic nature of the problem, but one that always appears in Fourier analysis, namely 
how well can one approximate a function if one uses only finitely many of its Fourier coefficients. 



2. Continuous functions on aperiodic point sets 

2.1. Local hulls. 

We work in R d , a real Euclidean space of finite dimension d. For r > 0, let D r (R d ), be all point sets 
A C R d for which the distance \x — y\ > r for all x, y € A. This means that D r consists of all discrete point 
sets with minimal separation > r between their points. 

The local topology on D r (M. d ) can be intuitively introduced as follows. Two sets Ai and A2 of D r (R d ) 
are 'close' if, for some large R and some small e, one has 

Ai n B R c A 2 + £? e , 

A 2 nB fl c A-i + B e , (1) 



More generally a compact Abelian group 
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where Br (resp. B e ) is the ball of radius R (resp. e) around 0. Thus for each point of Ai within the ball 
Br, there is a point of A2 within the distance e of that point; and vice versa. Pairs (Ai, A2) satisfying (p} 
are called (R, e)-close. 




B R (0) 



» points of A 1 
. points of A 2 

Figure 1. Two sets Ai and A 2 of D r (R d ) that are close. 
The local topology is actually a metric topology, although we make no use of this fact here. 
Definition 1. For A e D r (R d ) the local hull of A is 



X(A) = {t + A : t € R d } C D r (M. d ) , (2) 
i.e. take all translates of A and take their closure in the local topology. 

Proposition 2.1. [21 The translation action ofM. d on A lifts to a translation action on X(A). The local 

hull X(A) is compact and the M. d -action on it is continuous. □ 

Example 1. Let A be the lattice Z d in R d . Then X(A) = R d /Z d is the d-torus (with its usual topology). 
Intuitively we translate the lattice around with M. d . But translation by an element of Z d leaves A invariant. 
So R d /Z d parameterizes all distinct positions of A under translation. 

Example 2. Let A be any Penrose tiling. Then X(A) is the set of all Penrose tilings that are locally 
indistinguishable from some translate of A. X(A) contains considerably more than just the translations of 
A. In fact X(A) consists of all Penrose tilings based on the same pair of Penrose rhombs and the same 
orientations as appearing in A. 

Generally one may think of X(A) as some sort of local indistinguishability class of A. 

2.2. Continuous functions on A (A). We assume that A C D r (R d ) and AT (A) are as in W2.1{ 

Consider a function 

F:X(A) — ► C. (3) 

We can define from it a function 

f :R d — > C. 

by 

f(t)=F(t + A). (4) 
If F is continuous then we note that for all ti,t 2 € K d , 

ti + A, t,2 + A are close 

=^ F(t 1 + A), F(t 2 + A) are close 
fih) are close. 

Thus continuity of F implies continuity of /, and we see that / is local, or almost periodic, with respect to 
A in the following sense: 
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Definition 2. A function / : R d — ► C is called local with respect to a set A S Z? r (]R rf ), or A-local, if for 
all e' > there exist R and e so that whenever ti,t2 € R d satisfy that ti + A and + A are (i?, e)-close then 

|/( fl )-/(t 2 )|<e'. 

The intuitive meaning of this is that / has the very natural property, at least from the perspective of 
physical systems, that it looks very much the same at places where the local environment looks the same. 
Local functions are easily seen to be continuous on R d . 

Using locality, we can go in the opposite direction. Let A £ D r (R d ) and let / : R d — > C be local with 
respect to A. Define 

F : {t + A : t € R d } — ► C 
(so that F is a function on a part of D r (R d )) by 

+ A) = /(*). 

Then F is continuous on {t + A : f e R" 1 } with respect to the local topology. In fact it is uniformly 
continuous. The reason for this is that the continuity condition which defines the local-ness of / is based on 
the uniformity (i.e. the notion of (R, e)-closeness) defining the local topology of {t + A : t E R d }. 

It follows that F lifts uniquely to a continuous function ([3]) on the local hull A (A). 

Proposition 2.2. For each local function f with respect to A there is a unique continuous function ([3]) on 
the local hull, whose restriction to the orbit of A is f . Every continuous function on the local hull of A arises 
in this way. □ 

Thus we see that a locality with respect to A and the existence and continuity of an extension function 
on X(A) amount to the same thing. 

In the situation that A(A) is equipped with an R d -invariant probability measure /i (i.e. a positive Borel 
measure \i with /x(t + A) = n(A) for all Borel sets and with fi(X(A)) = 1), the action of R d on X(A) leads 
to unitary action T of R d on L 2 {X(A),fi). Namely for all F € L 2 (X(A),/i) and for all t € R d , T t F is the 
function defined by T t F(T) = F(-t + T) and with 



(F | G) := [ FGdfj, 
Jx(A) 



A' (A) 

we have 

(T t F | T t G) =(F\G). 

In principle the spectral theory of L 2 (X(A), /x) should allow one analyze A-local functions / : R d — » C by 
analyzing their corresponding functions F on L 2 (X(A), /x). 

For one very important class of subsets A this can actually be carried out in detail - namely the class of 
model sets, which we now introduce. 

3. Local functions on model sets 

3.1. Cut and project schemes, model sets and torus parametrization. An important class of point 
sets A C R d for which we know a considerable amount about the corresponding hulls A (A) is the class of 
cut and project sets, or the model sets as they are often called [TBlflT] . 
Consider the cut and project scheme 

R d J\_ R d x R d R d 

U (5) 

1 — 1 ~ dense image 

and a window Q. Here L is a lattice in R d x R d which is oriented so that the projections into R are 1 — 1 
and dense respectively. 

In ([5|) the left-hand R d is physical space, the space in which A is going to lie. The right-hand R d is 
internal space, the one that will be used to control the projection of the lattice L into physical space. The 
image of L under projection into physical space is denoted by L. Since this projection is one-one, L ~ L as 
groups, so L is a free Abelian group of rank 2d, i.e. it has a Z-basis of 2c? elements. However, it necessarily 
has accumulation points, and the typical situation is that L is dense in physical space. 
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It is convenient to use notation like x, x, x' for the elements of L and their respective left and right 
projections. Then x = (x, x') where x runs through L. This implies existence of the mapping ( • )' : L — > L' 
defined by x i— > x' , which passes from physical to internal space. 

Note that this mapping ( • )', as given, is only defined on L. It cannot be extended in any canonical way 
to a mapping K d -> l d . However it does extend canonically to the rational spans of the objects in question, 
and we shall make use of this later. 

We choose a subset ft in internal space. This window is assumed to be compact, equal to the closure of 
its interior, and to have boundary of measure 0. Using it we define 

A = A(Q) = {x : x G L, x' G fi} . (6) 

Sets of the form t + A(f2), t G R d , are called cut and project sets or model se£s0. In particular, for each 
(x, y) € R d x R d we may define 

A( B|W ) = x + A(-y + fi) . (7) 

If (x, y) = (x' , y') mod L, then x + A(—y + fi) = x' + A(—y' + f2), as can be verified directly from the 
definitions. Thus these model sets are parameterized by the torus of dimension 2d, 

(R d x R d )/L ~ (R/Z) d = : T . (8) 

We note, though, that the parametrization need not necessarily be 1 — 1, i.e. in general, A-( x ,y) — ^(x'.y 1 ) 
(x,y) = (x',y') mod L. In the sequel, for simplicity, we shall usually write (x,y)z, for the congruence class 
(x, y) mod L. 

There is a natural measure, the Haar measure, 9j on T. This measure is the obvious 'area' measure in 
the case of T 2 , and 'length' measure for T 1 . It is invariant under the M d -action: R d acts on R 2d /L by 

t+(x,y) L = (t + x,y) L . (9) 

Of particular importance to us is natural embedding (see Fig. 2.) 

R d — > T t»(t,0) L , (10) 

which lies behind the connection between almost periodicity in physical space and periodicity in some higher 
dimensional setting. The image of this mapping is easily established to be dense in T. 

Now start with A = Am,o) = A(f2), and translate it by elements t G R d : 

t + A(fi) = t + A(0 + O) = A (t)0) . (11) 
Form the local hull X(A), the closure of the set of all translates A( t ) of A under the local topology ([2|). 
Proposition 3.1. |22j There is a continuous mapping, 

P:X(A) — > T, (12) 

called the torus parametrization, such that 

1) P is onto; 

2) fi is 1 - 1 almost everywhere, seen from the perspective of Haar measure on T; 

3) for all t e R d , for all A' G X{A), one has /3(t + A') = t + j3(A'). 

4) f3(t + A) = (t, 0) L for allteR d . □ 

X(A) and T are both compact spaces with natural R d action, and /3 : X(A) — > T is an onto R d mapping. 
But X(A) and T are subtly different. Although t + A G X(A) and t + A i— > Au m for all t, not every 
element of X(A) is a A-( x ,y) for some (x,y). Rather, when [—y + <9f2) n L' ^ 0, then for all x G R d there 
are always Ai ^ A2 of elements of X(A) that are mapped by (3 to the same (x, y)L- When z = (z, z') G L 
with z' G (-y + then there will always be Ai, A 2 G X(A) with /3(X(Ai)) = /3(X(A 2 )) = (x,y) L , yet 
x + z e Ai, x + z ^ A 2 . In other words, there are ambiguities regarding the lattice points that project onto 
the boundary <9f2 of il. An example of this is shown in the footnote appearing in §5.31 The fact that (3 is 
1 — 1 almost everywhere is due to our assumption that the boundary of f2 has measure 0. 

When we say that (3 is 1 - 1 almost everywhere, we mean that the set A of points £ in T, for which there 
is more than one point set over £, satisfies 9j(A) = 0. 

There is a unique R d -invariant ergodic measure /i on X(A) with /x(X(A)) = 1. In fact (3 relates fi and 9: 

f3(fi) = d. (13) 



2 Model sets can be taken more generally with any locally compact Abelian group as the internal space. 
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Or more specifically, 6(A) = for all measurable subsets A of T. With in hand, we can introduce 

the space L 2 (X (A) , /j,) of square integrable functions on X(A). As already pointed out, the natural action 
of M. d on this is unitary. 

3.2. From hulls to tori. 

Square integrable functions on X(A) and square integrable functions on T can be identified, 

L 2 (X(A)^)^L 2 (T,6). (14) 

The isomorphism is easy to understand: 

X(A)^-^J 

\ 

x G 
Go/3 N 

AT 

C 

gives us a map 

L 2 (T,8) — » L 2 (X(X),ii). 
Since (3 is almost everywhere 1-1, the map is a bijection. 

This allows us to analyze functions on X(A) by treating them as functions on T. The advantage of this 
is that functions in L 2 (T, 9) have Fourier expansions 

F(z) = a k e 2wi{%lz) , zeT, a k eC. (15) 

Here (•, •) is the natural dot product on R d x M d = M. 2d and L° is the lattice which is Z-dual to L: 

keL° (k | x) € Z , for all x £ L . 

Assuming that the dot product is rational valued on L, this new lattice is in the Q span of L and its 
elements k have the same type of decompositions k = (k, k') as elements of L. In particular, for each k there 
is a unique k £ R d . We prepare for the ultimate reduction of everything to the physical space by already 
using the symbols rather than for the coefficients of the Fourier expansion. We shall write L° for the 

projection of L° on the physical side, so we can just as well write k S L° as k £ L° . We will call L° the dual 
(or Z[r]-dual) of L. The mapping (•)' : L° — > R d here is compatible with the one on L, indeed QL° = QL. 

The corresponding functions on X(A) have similar expansions. This works as follows: If F : X(A) — > C 
corresponds to F : T — > C then, for Ai S A(A) with /3(Ax) — (x, y)^, we have as L 2 -functions 

F(A X ) = F((x, y) L ) =Y, a * e 2 ^ % \^ . (16) 

k£L° 

3.3. Local functions on model sets. 

Let A be a model set arising from the cut and project scheme of (0. Let T = T 2d = (R d x K d )/L be the 
torus with torus parametrization 

(3:X(A) — ► T. 

Then we have the identification ([M]) of the corresponding L 2 spaces. Each element A' e A (A) maps by f3 
to a point /3 (A') in T. We also know (3(A) — (0, 0)^, and (3(t + A) — > (t, 0)^. So we know how /3 works on 
R d + A. 

Suppose / is a local function with respect to the model set A. From the local function / we have its 
extension F £ L 2 (A(A),/i) which is continuous. Then we obtain F £ L 2 (T,8), where 

F((t, 0) L ) = F(f3(t + A)) = F(t + A) = /(*) , (17) 

and we can write 

F(-)= J2a k e 2 ^, 

where L° is the dual of L. 

Definition 3. The Fourier-Bohr expansion of the local function / is 

f(t) = F(t + A) = F((t, 0) L ) = ak e 2 ™ <S|(t ' 0)> = ak ■ ( 18 ) 
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Figure 2. A fragment of the orbit of A as seen in the torus parametrization. 



Our study of almost periodic functions on A becomes the study of functions on T and their restrictions 
to the 'spiral' orbit R d + A in T given by the embedding (fT0|) . i.e. restriction to (R d , Q)l- 

3.4. Fourier coefficients. 

Let us continue with the situation in §3.31 For all x € T, 

F(x) = afce 27 "^ (19) 

k£L° 

where a k = [ e~ 2m{ ^ x) F(x) d9 T {x) . (20) 

Unfortunately, we don't have total control over F. We know it only on (M' 1 ,0)l. To compute au out of / 
alone, we use the Birkhoff ergodic theorerr|| : for all continuous functions G on T, 

G{x)d6 T {x) = lim — l — f G((t,0) L )dt. (21) 
R-kx> voliJR J Br 

Thus 

a k = hm — L- / e- 2 ^~ k ^F((t,0) L )dt 

fl->oo vol B R J Br 

= lim -J— [ e- 2 ^f(t)dt. (22) 

R->oo VOi Br J Br 

Here we use F((t,0) L ) = f(t) and k — (k,k'), so 

(k\ (t,0)) = (k,t) + (k',0) = {k,t). 

The averaging sequence {Br} that we have employed here can be replaced by any unbounded ascending 
sequence {A„}, where A n C M. d is compact, 1J A n = R d , and the boundary of A n has measure for all n. In 
actual practice one should adapt the averaging sequence to the problem at hand. 

4. Discretization 

4.1. Main components. Our objective is to devise a discrete method by which to estimate the coefficients 
of (|22| of a local function / with respect to a model set A = A(f2). This means replacing the integral 

by a finite sum of values of the integrand. The problem is to do this in such a way that it respects the cut 

and project scheme in which the model set lives, can be guaranteed to converge in the limit to the required 

integral, and can be carried out efficiently from a computational point of view. 
There are three components to this: 

(i) deciding on a suitable domain of integration (what should we use for Br?)] 

(ii) creating the points of evaluation of the integrand, including how many there should be; 

(iii) deciding to which set of Fourier coefficients (which values of k) we should restrict our attention. 



^It is hard to find convenient references for the Birkhoff ergodic theorem in the form we need it. Most references prove it 
over Z or R. In 1 1 21 there is a proof over Z d which is easy to adapt to M d . 
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A key feature of discrete methods involving periodic functions is the use of finite groups arising from 
refinements of the period lattice and quotients of its dual lattice, e.g. [20]. We need to translate this concept 
into the context of cut and project schemes. The set of points on which the integrand is evaluated is created 
out of the same cut and project process that creates the original model set. The ingredients are a choice of 
a suitable lattice Ln D L, which then gives rise to the finite group Ln/L. The data points in M. d at which 
computations of our functions will be made come by projection into physical space of a suitable set of coset 
representatives of Ln modulo L. The corresponding frequencies (wave vectors) k are chosen from the dual 
lattice L° . The choice of values of k at which we should evaluate the Fourier coefficients come by selecting 

— o 

suitable representatives of L° modulo L jv • The key point is the duality 

1 

N~ 



:L°/L° N xL N /L^ -Z/Z. 



4.2. Outline of the discretization process. In this section we give more precise details as to how the 
goals of i j4.ll can be achieved. The process necessarily involves a number of decisions, which can only be 
made in the context of the situation at hand. In $5]and $5] we shall see how this looks in particular examples. 

It should be noted that although getting the computational details set up is somewhat involved, these 
details depend only on the cut and project scheme and the degree of accuracy required from the computation. 
Once this is established the data points and choices of frequencies are already determined, and they suffice for 
the Fourier analysis of all functions that arise out of the same almost periodic family, and are algorithmically 
easy to compute. 

We begin with the cut and project scheme ([5]) with torus T and note the natural extension of the mapping 
(•)' to the rational span of the module L: 

QL QL — ► QL' (23) 

X < — I x = (x, X ) I— > X 

We assume that R 2d ~ M. d x M. d is supplied with the standard dot product (denoted (-|-)), and then define 
the dual lattice: 

Z° = {Y € R d x R d : (Y | x) € Z for all x € L}. 

Then L° is a Z-module of the same rank as L, namely 2d. There is a cut and project scheme dual to §5§ 
of which L° is the lattice H71: 



— m d x 



jj> 1-1 Jo den se im age (L°\' 



(24) 



We shall use the same type of notation as in (j2li|) for this scheme too. It arises by taking the Pontryagin 
duals of all the groups in (|2"5|) , whereupon L° appears as the dual of the torus T. For more on dual cut and 
project schemes see [17] . 

Dualizing can be considerably simplified if the inner product on R d x R d is rational valued on the lattice 
L. This quite often happens in actual practice. For instance, it happens in the Fibonacci example below, 
where the inner product arises from the trace form on Q[\/5]- When this happens one can identify L as a 
subset of QL and thereby avoid having to find the new (•)' mapping. However, in the general situation such 
simplifications need not exist, and we do not assume them here. 

In reading what follows it is good to keep in mind the equation 

e 27ri(k\S) _ e 2ni(k\s) e 2Tri(k'\s') ^25) 

which lies at the bottom of the approximation. The discrete Fourier analysis is accomplished by a dual pair 
of finite groups from which the k and the s will come. The actual values of k and s are important only 
modulo the lattices Lm and L respectively, and this freedom lies at the heart of the process. 

The values of s should be in the range A of our integration, and ideally we would have the corresponding 
s' = since the s are supposed to be representing points of the physical space R d . However, the latter is 
not possible, so we attempt to choose the s G A along with s' as small as possible. The approximation then 
works by throwing away the term which involves k' , s' in (|25[) . The set of values of k is constrained primarily 
by the requirement that the values of \k\ should be small, see the discussion after Step 5 below. 
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4.3. The six steps. Step 1: Choose a finite subgroup of T of order N. This appears in the form Ln/L 
where Ln D L is a lattice refining L. The number N will determine the number of points of evaluation in 
approximating the integrals by sums. We have Ln C QL and it has a Z-dual L° N C L° which is of index N 
in L°. 

This affords the natural pairing 

L°/L° N xL N /L — > iz/Z (26) 

induced by (• | •) (and still denoted by (• | •)). Using the (•)' mappings, we obtain, in the obvious notation, 
the Z-modules Ln D L , L° N C L°, Ln/L ~ Ln/L, L° / L° N ~ L/Ln, and an induced pairing 

(• I •> : L°/L° N xL N /L lz/Z. 

Replacing L by Ljy, we have the refined cut and project scheme 

R d < — R d x R d — > R d 

U (27) 

ijy < > £jV < ► ijy 

and similarly its dual. 

Step 2: Choose a fundamental domain C for L. A canonical choice would be the Voronoi cell of L at 0, 
but any other choice is allowable. In the examples below we use the parallelogram defined by a pair of basis 
vectors of L. Cover (R d ,0) C R 2d with a set of translates t + C of C by elements of L. The projection of 
these cells into the physical space R d covers it, though in general the projected t + C have many overlaps. 

Step 3: Form S := Ln H C. This provides a complete set of representatives in Ln for the group Ln/L. 

Step 4: Choose a region A which will be the delimiting the range over which the Fourier coefficients of / 
will be estimated in the form 

a k = ^ i/;-^ /(*)*. 

These integrals have to be computed for values of k that come from L° /L° N . 

Take as A the image of a finite set of the translates of C appearing in Step 2, i.e. 

A = {{Jt + C)W 

teT 

for some finite subset of T elements of L. The choice of A is again determined by the problem at hand. 

We need next to determine a set of data points in R d which will serve to replace the integrals of Step 4 
by finite sums. This is the purpose of the next step. 

As we pointed out above, we are free to translate the elements of S by L as we please, and we wish to do 
this so that the projected images arc in our region of integration A = [j^ e f(t + C)". We also keep in mind 
that we wish to do this so as to minimize the size of the corresponding s'. 

Step 5: For each s — (s, s') <G S find a t(s) G T for which \s' +t(s)'\ is minimal. Then the set of data points 
is 

D := {s + t(s) :seS}. 

At this point, for each k e L° we have 

ueD 

and the resulting approximation of / is 

keK 

The set K is to run over a complete set of representatives of L° /L° N . The choice seems free, but one may 
assume that in most cases the lower frequencies (smaller |fc|) are most essential in approximating / by using 
only finitely many of its frequencies. For this reason we have: 

Step 6: Choose k for each class of L° /L° N with |fc| taken as small as possible. 
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This concludes the broad description of the algorithm. 

5. A Fibonacci example 

To make all this more concrete we work through the details of a one dimensional example, the well-known 
Fibonacci sequence, where the cut and project scheme lives in two dimensions and the geometry of the data 
set D and the set of translates T via a refinement lattice and new windows are easily visualized. This involves 
first setting up the cut and project scheme in detail, £ j5 . 1 1 and £15.21 and then describing a Fibonacci point 
set arising from the standard Fibonacci substitution in terms of it, ^5.31 We then follow Steps 1 through 6 
of §4.3[ which provide the data points and corresponding frequencies that will work for the analysis of any 
local function that we may choose. In $6]we apply this information to two simple examples of local functions 
to see how well the methods actually work. 

One of the great virtues of the cut and project method is that it is primarily an algebraic tool and does not 
require great geometric insight to use it. Given that for aperiodic structures in dimension greater than 1 we 
almost always in the situation of lattices of rank greater than 3, and hence forced into spaces in dimensions 
greater than 3, this type of algebraic formalism is of enormous value. However, in the Fibonacci example, in 
which everything can be done in 2-dimensions, it is useful to see the underlying geometry explicitly. Thus 
we have gone to some effort to show the geometric meaning of the central feature of the method, that is, 
the creation of the data points, and to show what the approximations look like and how good they are in 
comparison with exact computation. In actual practice all this is unnecessary. The only part of the algorithm 
that requires any serious insight into the geometry (and it is actually trivial in the Fibonacci example) is 
the selection of the fundamental domain and the translates of it that are to be used. It is their projection 
that makes the domain in physical space in which the data will lie. 

5.1. The Fibonacci cut and project scheme. Let t := i(l + ^/5) and let Z[t] = Z + Zr. Then Z[t] is 
the ring of integers of the field Q[r] = Q[\/5] and r 2 = r + 1. Let (•)' on Z[t] and Q[r] be the conjugation 
that interchanges and — V5- 
We define 

Z[r] := {(x,x') : x £ Z[t]} clxl. 
Z[t] is a lattice inlxl and its natural projections 

into K provide the set-up for the Fibonacci cut-and-project scheme: 

M < — MxM — ► R 

u u u 

Z[t] <— zjrj — ► Z[t] (28) 

x < — i x = (x,x) I ^ x' . 

A natural basis of Z[r] is {(1, 1), (t, t')} and the standard inner product on Z[t] is defined by 

(x\y) = (2x- y) z = Sxy + (Sxy)' , (29) 

Here the notation ( • )z indicates taking the rational component a of 2x ■ y € Z[t] =: a + br, and 

«y=(r>/5)- 1 = (r(r-^))- 1 = ^ I) 

5' = (-r'VS)- 1 = (t'(t' - r))- 1 = p^TT > ( 3 °) 

In particular, 

(T|T) = (2(1-1)) Z = 2, (r |?H(2r 2 ) z = 2, (T | r) = (2r) z = . 

The geometry of the fundamental cell for the lattice Z[r] is illustrated in Fig. 3. More details of the material 
here can be found in [ID] , 
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5.2. The dual lattice. The inner product (• | •) allows us to identify Z[r] inside the rational span of Z[r] 
and to express (• | •) in terms of (■ | ■). The basis dual to {(1, 1), (r, r')} is given by 

u>i = \l = {\,\), £52 = £r= (§,£), (31) 

and the dual lattice is 

(Z[r])° = ZSi + Z£ 2 = ■ ( 32 ) 
The elements of (Z[r])° are always of the form k = (k, k') and we can write 

Z[T]°:={fc : k= (k,k') £ (zjr]) } (33) 

Then 

(Z[r])° = Z[t]° = : fc = (ft, fe') € Z[t]°} , 

and (ft | x) becomes (k \ x), which is more useful notation for the sequel. 
The 2-torus T of the cut-and-project scheme (|2"51) is then 

T=(Rx K)/Z[r] . 
Fourier series on T are expressed in terms of the characters x%i namely 

x ~(x) = e 2 ™W £ ) , where k = (k, k') £ Z2i + Z£ 2 . (34) 

Here x can be arbitrary in R x K, but we shall need to compute only with x £ QcDi + Qa;2 for which ~ is 
well-defined (|28|) . For these elements (fc | x) = (2k ■ x)q. 

The product (• | •) extends to R x R, but care has to be taken. For a typical element al + br, a, b £ K of 
our superspace IxK the inner product is calculated as 

((al + br) | (cl + dr)) = 2ac + 2bd , 

where 

al + br = a(l, 1) + b(r, t') = (a + br, a + br') elxl. 

Consider 

fc = aii + 6S 2 = fT+|?= (^i^r 1 -) =(k,k')£Z\r~}° , 

where a,b £ Z. Suppose we want to compute (k \ (t, 0)), t £ R. From 

( 1 '°)-775 T +7I ?: = « + ^' 
(0, 1) = t 2 S1 -t5t =6'1 + t'S't , 

where S, 5' are from (|3U)l . we have 

(k | (1,0)) = 25k, 

and so, 

(k | (t,0)) = 26kt, for all (€». 

We are interested in model sets A coming from the cut and project scheme (f2"5|) . In the context of (|2"5)) . the 
Fourier-Bohr expansion (|18[) assumes the simpler form 

/(f)=F(t + A)=F((t,0)^)= J] « fe e 2 - 25fct . (35) 

fcGZ[-r]° 

We shall use the au computed in the form 

a k = lim I [* e-™ 2kt f(t)dt. 



For future use note that 



(fc | (0, 1)) = ((2^, 2^) | S'l + t'S't) = 2k' 5' 
(k | (0,u)) = 2uk'S' . (36) 
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5.3. The Fibonacci model set. 

The standard 2-letter Fibonacci sequence is the fixed point of the substitution a — > a6, 6 — > a: namely, 

abaababa .... 

With tile lengths r for a symbols and 1 for b symbols, and starting at 0, we get the sequence of tiles that 
cover the non-negative part of the real line. The left-hand ends of these tiles, 

0, t, t + 1, 2r + 1, 3t + 1, 3t + 2, 4r + 2, 4r + 3, 5t + 3, . . . 

form an infinite sequence of points on the non-negative real line. This set appears explicitly as the non- 
negative part of the model set 

A = A([-1,±)) 

arising from the cut and project scheme ([28|) F1 A decomposes as A a = A([— ^)) and Ab = A([— 1, —75)), 
which give the left-hand end points of the a and b points respectively. 



Internal 
Space 


(0,1) 


1=0,1) 




±< 




/ b a v. 




1 J 

X 

Physical 




a b / 
(1,0) / 




space , 

f ~* 

T 




1 -f 

a 1 b / 


> 


\< 
% 




a | a / 








t=(t,t') 





Figure 3. The fundamental cell C of the lattice Z[r] = Z(l, 1) + Z(r, r') is shown in terms 
of the standard coordinate system for 1x1. When C is imagined as a torus by identifying 
opposite sides and the physical space is wrapped onto it using these identifications, then the 
a and b regions shown here indicate which parts of the physical space are lying in a or b 
tiles respectively. 

A useful and commonly used way to visualize the distribution of o and b points is to view them on the 
torus T after the embedding (fiT))) of M in T. The fundamental cell, see Fig. 3, with opposite edges identified 
is T. To see A itself being formed, we start at (0,0) and trace out (t, 0), ( £ R, with the usual rules for 
exiting and re-entering the fundamental cell. We move continuously. Our moving point has three 'whiskers' 
attached to it. 

a(l): a whisker of length i facing vertically upwards; 
a(2): a whisker of length facing down; 

6(1): a whisker of length i that faces down, but has an initial gap of size \. 

The rule is this. As our point moves along the line on the fundamental cell, if the a(l) whisker hits (1, 1) 
(i.e. (1, 1) is close enough to get cut by this whisker), we get an a-point of A. If the a(2) whisker hits the 
point (r, t'), we also get an a-point of A. If the 6(1) whisker hits the point (r, r'), we get a 6-point. 

Alternatively one can use A((— 1, —]), which differs from A in the two points coming from the ends of the interval [— 1, — ). 
This is an example of a pair of sets that map to the same place in T. This ambiguity shows up in an interesting way later on, 
see 323] 
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Essentially the window, represented by the whiskers, is carried along with the moving point, and exits 
and re-enters the fundamental cell with the moving point. 

The advantage of this point of view is that it allows us to divide the fundamental cell into regions which 
deliver the a- and 6-points, and thus to see what our function f(t) looks like on the fundamental cell. 

In Fig. 3 the vertical dotted lines are the key thing. Our moving point moves to the right and every 
time it crosses a dotted vertical line, we get a new point of A, thus starting a new interval. We stay on that 
interval until the next crossing. 

The a and b regions are indicated in the Fig. 3. 

5.4. Discretization. 

We now go into the details of how to deal with the Fourier analysis using discrete methods. Here we 
follow the six steps outlined in ^4.31 

Write L for Z[r], so L = Zl + Zr. Let C := {ul + vt : < u, v > 1} be a fundamental region for L. Its 
volume is equal to v5 . 

Fix any TV € Z + . We want to create a lattice L^ that is a refinement of L. In the £)4.3I we took Ljy D L 
with [Ln : L] = N. In our present situation we use the most obvious lattices, namely -wL. These actually 
have index N 2 , so there are some slight notational differences between this section and £)4.3l Since iVXjv Q L, 
so Ln C Ql + Qt. 

Let S :— Ln n C, so S is a complete set of representatives of Ljv mod L. 
Together this completes steps 1, 2, 3 of £ 14.21 

For understanding the approximations better we also introduce a fundamental region CV for Ln, chosen 



For Step 4 we choose as our delimiting range in R an interval [0, R] where R is a positive real number 
taken so that (R, 0) taken modulo L is on the boundary of C. Following along the path (t, 0), t S M>o, and 
wrapping around C as indicated in Fig. 4, this amounts to stopping at some point R where the path is just 
exiting C, so that we have an exact number of passes of C. Thus the path (t, 0), < t < R involves an 
explicit set of translates U + C, i = 1, . . . , M of C. Let T = {t%, . . . , 4m}- 

Let / : R — ► C be any continuous local function with respect to the model set A = A([— 1, ^)), and let 
F : T — > C be its extension to a continuous function on T = R 2 /L. Generally we are interested in Fourier 
decomposition of F and along with it the corresponding decomposition of /. Thus we wish to compute 
things like 



so that 




(37) 




Here _F(-)e -27 ™( fc l') is just some other continuous function on T that is local with respect to A. So it suffices 
to deal with some general continuous function G on T and its restriction g(t) = G((t, 0)) to the line (R, 0) 
mod L. 



Let 



en '■= sup sup \G{x) — G(si)\ . 

i=l,...,N 2 x£si+C N 



Our first estimate is 



Jt 

^f, the error in this is estimated 




N 2 



Since vol Cat 



by 




2 



2 



(38) 




14 



R.V. MOODY, M. NESTERENKO, AND J. PATERA 



i.e. the error in our approximation is bounded by v5£iV- In practice Ln and N would be chosen so as to 
provide a suitably a priori assigned value of en- 

5.5. Restriction to /. 

In order to be useful, the computation must be restricted to values of the function /, since this is all that 
one is given in practice. We wish to use an approximation of the form 

i [ g(t)dt = ^ [ R G(t,0)dt ~ [ddSr, 

ft JO ft J0 JT 

in accordance with (p?Tj) . 

In this section we indicate the geometry behind Step 5. Again, it should be pointed out that in the final 
analysis much of the detail that appears here need not appear at all in the actual algorithm. 
The path {(t, 0)l | < t < R}, wrapped around C, is divided by C into segments 1%, . . ., Im- 



Internal 
Space 



start (0,0), 



end (R,0) L 




Figure 4. Ten crossings of C by the path {(t, 0)i | < t < R} are shown in the torus. 
(Opposite sides of the parallelogram coincide.) 



Let bi < &2 < • • • < om be the projections onto internal space of the boundary cutting points of 
{(t,0) : < t < R}. Define 

J „ fci+&2 „ b M -l+bM i 

co = t , ci = 2 , • • • , cm-i = 2 ' CM = 1 • 

The strips Si formed by passing the interval [ci_i,c;) in internal space through C in the direction of the 
physical axis, form a partition of C. Through each strip in C runs a part of our line {(t, 0)l ■ < t < R}. 
The idea is to use the intervals [cj_i, Ci) as windows and the line segments U as (part of) the physical space 
for a model set construction based on the lattice L^. This will then produce the points on which will be 
our data points for the evaluation of the functions G, and then g. In other words, we are implicitly using 
the partial model sets Aj~(Si) n U- 
Let m(i?) :— max{|c; — c^_i | } and 



z N := sup sup 

i (x,v),(x,u)eSi 



G{x, v) — G{x, u) 



1 M. 



Each §j € S lies in exactly one strip Si. Let (pj,qj) be its projection onto the line segment k. Then 
(Pj,qj) = (uj,0)i for some Uj € [0,i?]. Thus we obtain {(ui,Q), . . . , (u N 2,0)} on our path {(i,0)i : < t < 
R}. 

We have the estimate (l38l). Moreover, 



i— N i— N i— N 

i—1 i—l i—1 

since 

G( Pi ,q i ) = G((u i ,0))=g(u i ). 
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l=c 




Figure 5. The projections 61,... ,6m into internal space of the boundary cutting points 
of the path. The points Co, . . . , cm mark the boundaries of the smaller windows into which 
the original window [r', 1) is partitioned. 




lattice points 
projected lattice points 



(T,T') 



Figure 6. The points Sj lying in the strip Si formed by the window \c%-\, c<) are projected 
onto the segment I4 of the path {(t, 0)l ■ < t < R} producing data points Uj, shown here 
as small crosses. 



Thus 



%/5 n2 VE n2 



V5 



since Sj and (pi, qi) both have second component in the same interval [cj_i, Cj). 
Combining (|38|) and (|39|) , we have 



< V5 (ejv + e'jy) • 
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This provides a method of estimating the integral J r GdO using only g on [0, R] along with well chosen 
points in the interval. The two parameters N and R control the estimates. 

In spite of the apparent complexity of strips, what is going on is easy to implement. For each Si = (s^, s^) 
there is a translation vector tj G T for which + t'-\ is minimal. The corresonding data point is then 
u t := Si + tj. This is Step 5 of 

6. TWO EXPLICIT EXAMPLES 

6.1. Computation of Fourier coefficients. In this section we apply the methods outlined above to two 
almost periodic functions / : R — ► M, both local with respect to the model set A = A([— 1, -)) of £15.31 
The first is the distance-to-the-nearest neighbour function 

f(t) = the distance of t to the nearest point of A . (41) 

This is a continuous and piecewise linear function that is local with respect to A. 
The second is the function 

I 1, if a; is in a long interval , . 

j : j(x) = < . . . (42) 

I — 1, if a; is in a short interval , 

which local, but only piecewise continuous with breaks wherever x switches from a long to a short interval. 

Graphs of these functions are shown (solid lines) in Fig. 8 and Fig. 13 respectively. 

Our objective is to see how well the approximations we have discussed compare with the actual functions 
when the calculations are done with specific choices of data points. 

We approximate each of the two functions by a finite number of terms of its Fourier-Bohr expansion Q18p , 
which in our setting here reads: 

/(.*)= J2 «* e " Wl) . * eR - ( 43 ) 

For the two functions / chosen here, it is easy to determine the corresponding functions F on the torus 
T = (1 x 1)/L and hence to compute the Fourier-Bohr coefficients a>k exactly by (TlU)) . see £16.21 . For the 
nearest neighbour function these are given explicitly in (|50p . 

On the other hand, the approximation depends on the choice of the refinement lattice Ln, which shall 
here always be of the form (1/N)L = (1/N)Z[t] (so the index [Ljv : L] is A 2 ). This determines both the 
points that will eventually be projected into our data points and also the values of k which shall be included 
in approximating the sum (|43|) . namely a set K of k = {k,k') € Z[r] chosen as representatives of the group 
Z[r] /NZ[t] dual to the group ((l/N)Z[r])/L. Our choice is take the k with |fc| as small as possible, so 
K consists of elements k € Z[r], one for each congruence class of Z[r] modulo AZ[t], chosen so that \k\ is 
minimal in its class. We denote by f exact {x) the finite series, taken from ((43|) . approximating f(x), 

f exact ^ _ 0fce !hri(*W ; (44) 

keK 

Next we replace the coefficients in f exact (x) by their approximations 

1 r R 

4 ni :=^y f(x)e- 2ni ^dx (45) 

following from (f22|) . and denote the resulting function by f mt (x). We shall use various values of R, all of 
which correspond to a set of complete passes across the fundamental region, as illustrated in Fig. 5. 

The integrals (|45[) are to be estimated by reducing them to finite sums where the integrand is computed 
on the finite set of data arising as projections {tti, . . . , u^} of the Af 2 points in C, as explained in £14.31 

1 N 2 

The resulting approximation to f exact (x) is denoted f sum (x). This is the approximation that we have 
been working towards. Written out in full it reads: 

2 

27ri(k\uj) I 27ri(k\x) 




f-(x) = > { — > 7(« 3 -)e~*" l * |tyJ > ts^W* . (47) 
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It is computed out of data points in R, is a finite sum of exponential functions whose frequencies come from 
the Fourier module of /, and utilizes discrete groups arising out of the periodic setting of the underlying cut 
and project scheme. By its form, f sum is almost periodic and approximates the function F everywhere on 
the real line as it lies wrapped around the torus. 

Thus we have four functions: /, f exact ^ f mt , j surn ^ a ll defined for all x G R. The three approximation 
functions depend on the number M of passes of the real line going through the fundamental region C (see 
§ 15.41 and 15.51 for details), and on the number iV 2 of lattice points of Lm found in C. Each of them restricts 
the full summation of the Fourier-Bohr expansion to the same finite set K of frequencies. They differ due to 
the ways in which the Fourier-Bohr coefficients are obtained: they are exact in the first case (or at least as 
exact as real computation on computers can be), are derived from the integral approximation in the second, 
and come from the finite sum approximation to the integral in the third. The calculations and graphs shown 
below allow one to compare these functions. The exact coefficients and their approximations by integrals 
are sufficiently close as to make little difference to the graphs, so the figures are restricted to comparing 

J jint jsum 

By way of comparison, there is one further Fourier approximant, of the type that we are obtained by 
periodically extending / from its values on a fixed finite interval / of R to the entire space R. Thus we 
introduce 

f cos (x) = 2^ a k cos(-kirx), where at — - I f(x) cos(-kirx)dx . (48) 
fc=o " ^° 

Needless to say, f cos cannot be expected to have much relationship to / outside /, and it doesn't. In fact, 
jcos ^Qgg no ^ appear to be a good choice even on the limited domain /. 

6.2. Computing the exact Fourier-Bohr coefficients. It may be of interest to indicate how we com- 
puted the exact Fourier-Bohr coefficients. We can reorganize Fig. 3, translating the a and b regions into single 
blocks (see Fig. 7), a procedure that is familiar from the klotz construction see §7 often used in studying cut 
and project sets, [TB"] . 

It is then straightforward to understand the corresponding function F on the 2-dimensional torus, i.e. a 
genuinely periodic function for which f(t) = F(t, 0) mod Z[r]. 
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Figure 7. The a and b regions of Fig. 3 are reorganized here into two rectangles (heavy 
lines). The function F on T arising from the local function (|4Tj) is supported on these two 
rectangles. It is constant on vertical lines. Along the horizontal axis, on the b rectangle 
the function values start at 0, increase linearly to 1/2 at the mid point of the rectangle and 
then linearly decrease to again. Similarly on the a rectangle they increase to r/2 and then 
return to 0. 

Furthermore, it is easy to compute the 'exact' Fourier coefficients of F. Let k = auj\ + bu>2 € Z[r] , where 
a, b £ Z. Then 

(k | (x, y)) =x(k\ (1, 0)) + y(k\ (0, 1)) = 2xkd + 2yk'5' , 
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where 



a + br 



a + br' 



The Fourier coefficient for k is 

a k = [ F{w)e- 2 ^ w) dw= f F{u,v) e- 2m ^ uUv7) du dv . 
Jt Jo Jo 



(49) 



Writing 
we have 



ul + vt — (u + vt, u + vt) =: {x, y) , 



(k\ ul + vt) = 2kxS + 2k'y8' . 
Making the change from the variables u and v to x and y, and using the definition of F , we obtain 



1 f° f° (1 



i V2 



VeI-1 So (2 



1 

X+ 2 



T 

X ~2 



Note that the factor comes from x = u + vt, y = u + vt' . 



1 T 
1 t' 



e -i^(kx8+k'yS') dx dy _ 

-Vs. 



(50) 



6.3. Shadows of singularities. Inspection of the aperiodic approximations shows that they are remarkably 
faithful to the originals. But there is apparently a strange fuzziness in the approximation in the interval 
[— t 2 , 0]. The explanation for this is quite interesting and shows that the approximation method we use here 
is sensitive to quite subtle qualities of the aperiodicity. 

The particular Fibonacci set that we have used in our example is singular, that is to say, it is one at which 
the torus parametrization is not one-one. As pointed out in footnote of 15.31 A is but one of two distinct 
sets in A(A) that map to the same point on the torus T under the torus map. The other one differs from 
A only in that it contains {— r} but not { — 1}. The Fourier analysis takes place on T and treats these two 
sets equally. However the nearest neighbour functions of the two sets are different. This is a difference that 
it immaterial to F, which our function f sum is approximating, but which the original function F can see. 
What we see is the Fourier analysis hedging between the two scenarios. 

Further refinements of the lattice will never improve the situation on this interval. However, if we had 
used an element of A (A) at which the torus map were 1 — 1 (and this is the case with probabilistic certainty 
if one chooses randomly from A (A)) then this phenomenon would not occur and the approximation would 
be uniformly valid over the entire space. All one need do is to shift the window so that its end points do not 
lie in the set Z[r], e.g. \ + [—1, -]. 



7. Final Comments 

We should point out that there is considerable scope for adapting the scenario sketched out here. It 
is possible to arrange things so that they are related to tilings. One elegant tiling method is the klotz 
construction [T31[T3]. One begins with the decomposition into Voronoi cells of the lattice L in M. d x M. d , and 
along with it the corresponding dual cell decomposition into Delone cells. For each pair of P, Q consisting of 
intersecting d-dimensional faces P and Q from a Voronoi cell and a Delone cell respectively, form the klotz 
Pll x Q 1 - in R d x R d . The set of these klotze form a tiling of R d x R d . Furthermore the intersection of (R d , 0) 
with this tiling produces a tiling of physical space. (If one chooses instead to do the projections the other 
way around, one gets a different tiling). Several of these klotze can be combined to form a fundamental 
region for R d x R d , as we see in Fig. 7, and this type of choice should lead to computational methods that are 
adapted to the these tilings. For instance, the region A of integration might well be chosen as the union of 
a finite number of tiles. For more on determining Voronoi and Delone cells in the context of high symmetry, 
see [19]. 

The method advocated here is based on the idea of local functions, the extending of them into the context 
of compact Abelian groups, and the discretization of the resulting Fourier analysis by the use of finite groups. 



ALMOST PERIODIC EXPANSIONS 



19 



In the ID setting that was explored in detail here, the only symmetries involved arise from translational 
symmetry (which is at the base of the almost periodicity). In higher dimensions, especially those of interest 
to the quasicrystal community, decagonal, icosahedral, or other symmetries appear. In these cases there 
are a number of ways of utilizing the symmetry to greatly improve the efficiency of the computation, in the 
same spirit as [T5]. As we have pointed out, the preparations required for this depend largely on the cut 
and project scheme and not so much on the actual model set involved. Fortunately the cut and project 
schemes for these settings are essentially canonical and their strong algebraic nature makes this program 
quite feasible. 

Finally, although we have not spelled it out here, the way in which finite groups and their duals are used 
here makes the process amenable to the technique of the fast Fourier transform. Details will appear later. 
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Appendix A. Numerical and graphical data for two exact examples of section [6j 

I.E., THE DISTANCE-TO-THE-NEAREST NEIGHBOUR FUNCTION GIVEN IN (41) AND 
THE LOCAL FUNCTION DEFINED BY (42) 



Table 1. Comparison of the calculated coefficients at in the approximations j exact ^ f 
and f sum based on the parameters of Fig. 10. 



k 




"fc 


~sum 
"fc 


-± - ±T 
2 2 


-0.1065-0.03668i 


-0.1065-0.0371i 


-0.1086-0.0581i 


1 

2 


0.0243+0.0287i 


0.0236+0.0292i 


0.0269+0.0711i 


-\ + \r 


0.0026+0.0153i 


0.0035+0.0155i 


0.0233-0.0332i 


-¥ 


-0.0683+0.0407i 


-0.0680+0.0412i 


-0.0517+0.0542i 





0.3618 


0.3618 


0.3367 




-0.0683-0.0407i 


-0.0680-0.0412i 


-0.0517-0.0542i 


1 - At 

2 2 ' 


0.0026-0.0153i 


0.0035-0.0155i 


0.0233+0.0332i 


1 

2 


0.0243-0. 0287i 


0.0236-0.0292i 


0.0269-0.0711i 


2 ~ 2 


-0.1065+0.0367i 


-0.1065+0.0371i 


-0.1086+0.0581i 



Table 2. Comparison of the values of the function / of (|4"Tj) and its approximants /' 
f mt , f sum , and f cos , see Fig. 10. 



Xi 


/(*<) 


f exact {xi) 


f int (xi) 


f sum {xi) 


/ cos (^) 


-100 


0.8065 


0.6916 


0.6965 


0.7728 


0.1859 


-50 


0.4033 


0.4229 


0.4208 


0.4562 


0.3690 


-15 


0.3262 


0.2555 


0.2522 


0.1461 


0.4912 


-3-5t 





0.0577 


0.0584 


0.1378 


0.5365 








0.0658 


0.0670 


0.1165 


0.1859 


T 





0.0797 


0.0788 


0.0946 


0.1858 


0.25 + r 


0.2500 


0.2060 


0.2049 


0.1995 


0.3065 


0.5 + r 


0.5000 


0.3318 


0.3313 


0.3416 


0.3138 


1 + r 





0.1659 


0.1649 


0.1115 


0.3000 


1 + 1.25r 


0.4045 


0.3325 


0.3287 


0.1467 


0.5022 


1 + 2.5t 


0.8090 


0.6562 


0.6609 


0.6949 


0.4681 


1 + 2.75r 


0.4045 


0.3119 


0.3152 


0.2659 


0.2659 


50 


0.4033 


0.3265 


0.3229 


0.1209 


0.3690 


100 


0.1885 


0.3006 


0.3004 


0.2282 


0.1859 


500 


0.4396 


0.4669 


0.4651 


0.5364 


0.1859 
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FIGURE 8. The solid line is the graph of / of 1(21"]). The dotted curve shows the f sum 
approximation and the dashed curve shows the f mt approximation. Both approximations 
were calculated using a total of TV 2 = 49 lattice points in C and M = 11 path passes in C, 
corresponding to an interval of integration R « 23, 30. 
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Figure 9. The function / and approximants of (|41 [) are drawn with the same conventions as 
in Fig. 8, but for refined parameters: M = 17, TV 2 = 121, and R « 37, 43. See subsection [ITU 
for an explanation of the misfit of approximations in the region [— r 2 , 0]. 




Figure 10. The function / of (|4Tj) is graphed for parameters: TV 2 = 9 lattice points from 
C, M — 10 passes of the fundamental domain, and corresponding interval of integration 
given by R « 21,64. Also the periodic approximation f cos is shown. The number of 
cosine terms used is n = 50 and coefficients are calculated by formula (|4"8")) . Here the 
dotted curve corresponds to the f sum approximation, the dashed curve corresponds to f exact 
approximation, and the thin line represents f cos approximation. 
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Figure 11. Here we graph the same approximations as in Fig. 10, but the functions are 
now drawn on the interval [200, 215]. Observe how the aperiodic approximants continue to 
follow the graph of / while the periodic approximant now has no relation to it. 
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Figure 12. The approximations are presented with the same conventions as on Fig. 10 but 
the functions are drawn in the interval [—115, —100]. 

Table 3. Comparison of the calculated coefficients in the approximations f exact ^ / mt , 
and f sum . The parameters are as in Fig. 14. 



k 




a k 


~sum 
a k 


_I - I T 

2 2 


-0.1065-0.0367i 


-0.1065-0.0371i 


-0.1086-0.0581i 


_ 1 


0.0243+0.0287i 


0.0236+0.0292i 


0.0269+0.0711i 


-i + |r 


0.0026+0.0153i 


0.0035+0.0155i 


0.0233-0.0332i 


-b 


-0.0683+0.0407i 


-0.0680+0.0412i 


-0.0517+0.0542i 





0.3618 


0.3618 


0.3367 


1- 


-0.0683-0.0407i 


-0.0680-0.0412i 


-0.0517-0.0542i 


2 2 ' 


0.0026-0.0153i 


0.0035-0.0155i 


0.0233+0.0332i 


1 

2 


0.0243-0.0287i 


0.0236-0.0292i 


0.0269-0.0711i 


2 ' 2 ' 


-0.1065+0.0367i 


-0.1065+0.0371i 


-0.1086+0.0581i 



Table 4. Comparison of the values of the function / of (|42|) and its approximants /' 
f int , f sum , and f cos . 







f exact ^ 


f in \xi) 


f sum {xi) 


f cos (^) 


-100 


1 


1.0960 


1.0796 


1.1907 


0.8024 


-50 


1 


1.0690 


1.0628 


0.6269 


-0.1948 


-15 


1 


1.1020 


1.1637 


1.1570 


0.8162 


-3-5t 


1 


0.1092 


0.0963 


0.3515 


0.8093 





1 


0.5331 


0.5268 


1.5952 


0.8024 


T 


-1 


-0.0271 


-0.0291 


0.7757 


0.3036 


0.25 + t 


-1 


-1.2750 


-1.2957 


-1.3115 


-0.1820 


0.5 + r 


-1 


-0.7919 


-0.8123 


-0.7404 


-0.1884 


1+T 


1 


0.0704 


0.0971 


-0.0128 


0.7910 


1 + 1.25r 


1 


0.9261 


0.9625 


1.1215 


0.8005 


1 + 2.5t 


1 


1.1233 


1.1382 


1.0246 


0.7991 


1 + 2.75r 


1 


0.9577 


0.9519 


1.2470 


0.7826 


50 


1 


0.9241 


0.9374 


1.1735 


-0.1948 


100 


-1 


-1.2440 


-1.2652 


-1.5375 


0.8024 


500 


1 


1.0950 


1.0639 


0.9970 


0.8024 
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Figure 13. The solid line is the graph of / of (|42|) . The dotted and dashed curves 
correspond to the f sum and f mt approximations. Both approximations were calculated 
using a total of N 2 — 81 lattice points in C, M ~ 17 passes of C, and corresponding 
interval of integration given by R ~ 37, 01. 




FIGURE 14. The function / of (O is drawn for parameters: iV 2 = 49, M = 11, R w 23, 30. 
Also the periodic approximation f cos is shown as thin line. The number of cosine terms 
used is n = 50 and coefficients were calculated by (|48|). Here the dotted curve corresponds 
to the f sum approximation and the dashed curve corresponds to f exact approximation. 
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Figure 15. Here we adduce the same approximations as in Fig. 14, but the functions are 
drawn on the interval [200,215]. 
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Figure 16. The approximations are drawn with the same conventions as on Fig. 14, but 
on the interval [—115, —100]. 
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